Single grain heating due to inelastic cotunneling 
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We study heating effects of a singie metafile quantum dot weakly coupled to two leads. The dom- 
inant mechanism for heating at iow temperatures is due to inelastic electron cotunneling processes. 
We calculate the grain temperature profife as a function of grain parameters, bias voltage, and time 
and show that for nanoscale size grains the heating effects are pronounced and easily measurable in 
experiments. 
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Great efforts in contemporary materials science re- 
search focus on the understanding of thermal properties 
of nanoscale devices^^. The simplest possible nanode- 
vice consists of a single grain coupled to two leads (see 
Fig. [I} . The transport properties through this system are 
well understood*^ However, much less is known about 
heating effects in this system. The understanding of heat- 
ing effects on properties of nanodevices is especially im- 
portant for practical applications. Indeed, recent exper- 
imental researc h 8 ' 9 ' 10 has focused on the application of 
the thermal properties of low-dimensional devices for ef- 
ficient power generation. This defines an urgent quest for 
a quantitative description of heating effects on properties 
of nanodevices. 

Indeed, the effect of heating is very important for low- 
temperature nano-thcrmoelectric devices. In particular, 
next-generation devices made from nano-granular mate- 
rials can be influenced by internal heating, and it is im- 
portant to know the device's internal temperature. Gran- 
ular thermoelectric materials have the advantage that 
one can control the system parameters and therefore the 
device properties. A measure for the efficiency of a ther- 
moelectric material is the dimcnsionlcss figure of merit, 
ZT, where T is the temperature, which depends on the 
thermopower or Seebeck coefficient, S and the electric, 
a and thermal, K conductivities, ZT = S 2 aT/ kJ^ One 
has to take into account the fact that all three kinetic 
coefficients, a, K, and S, depend on the internal device 
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Recently we investigated the prop- 



erties of metallic granular metals in the metallic*^ and 
weakly coupled regimes^ which are relevant for appli- 
cations as low-tcmpcraturc thermocouples. Here we will 
answer the question to what extend heating effects influ- 
ence the properties of a single quantum dot - the building 
block of granular metals. 

In this paper we investigate the effect of inelastic co- 
tunneling on heating of a grain weakly coupled to two 
leads (see Fig. [1]). We calculate the grain temperature 
profile as a function of grain parameters, bias voltage, 
and time. For nanoscale size grains it is shown that the 
heating effects are pronounced and easily measurable in 
experiments. Here we consider the case of an isolated 
grain in the sense that it is not thermally coupled to an 



environment other than the leads, e.g. the grain is not 
coupled to a substrate. 

A single grain is characterized by two energy scales: 
(i) the mean energy level spacing S and (ii) the charging 
energy E c . We concentrate on the case of metallic grains 
which satisfy the inequality E c 3> 5. 

In addition to the above energy scales the system under 
consideration is characterized by the tunneling conduc- 
tance gt- In this paper we concentrate on the regime 
of Coulomb blockade corresponding to a weak coupling 
between the grain and the leads, gt <C 1. Our considera- 
tions are valid for temperatures S < T < E c . 
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FIG. 1: (Color online) Sketch of a single grain weakly coupled 
to two leads, g t <C 1. Temperatures T and T g refer to the 
leads and grain temperatures, respectively. Arrows indicate 
the inelastic cotunneling process through a single grain with 
a typical energy of the electron- hole pair T g . The dashed 
line corresponds to the Fermi level in the grain. Here only 
the electron propagation process responsible for the heating 
of the grain is shown (a- m ). The heat is removed from the 
grain by a similar process but with the electron starting and 
ending in the same lead (inelastic cotunneling loop, Ki n ) and 
by phonons (n p h). 



Transport mechanisms. — The only transport mech- 
anism contributing to grain heating is due to inelastic 
cotunneling. Cotunneling allows simultaneous charge 
transport through several junctions by means of coop- 
erative electron motion. Single cotunneling (Fig. [IJ, in- 
troduced in Ref. Q provides a conduction channel at low 



2 



applied biases, where otherwise the Coulomb blockade 
arising from electron-electron repulsion would suppress 
the current flow. The essence of a cotunncling process 
is that an electron tunnels via virtual states in interme- 
diate granules thus bypassing the huge Coulomb barrier. 
This can be visualized as coherent superposition of two 
events: tunneling of the electron into a granule and the 
simultaneous escape of another electron from the same 
granule. There are two distinct mechanisms of cotunncl- 
ing processes, elastic- and inelastic cotunneling. Elastic 
cotunneling means that the electron that leaves the dot 
has the same energy as the incoming one. In the event of 
inelastic cotunncling, the electron coming out of the dot 
has a different energy than the entering electron. This 
energy difference is absorbed by an electron-hole excita- 
tion in the dot, which is left behind in the course of the 
inelastic cotunncling process, Fig. [TJ Below we concen- 
trate on inelastic cotunncling because only this transport 
mechanism contributes to heating effects. In particular, 
elastic cotunneling and sequential tunneling do not create 
electron-hole pairs in the grain. 

In the following we consider the electric conductiv- 
ity due to inelastic processes, which heat the grain, and 
the electronic and phonon thermal coupling to the leads, 
which remove excess heat from the grain, for the descrip- 
tion of the complete heating effect of a single grain. 

Grain heating. — The grain heating due to inelastic co- 
tunncling can be described by the following kinetic equa- 
tion 
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where c v = (l/3)T g is is the heat capacitance of the grain 
with v being the density of states in the grain; dT g /dt 
describes the change of the grain temperature T g in time 
t; Cm is the cotunneling conductance (defined below), 
V being the bias voltage and a the grain size; and K is 
the thermal conductivity, with k = n ln + Kph where «i„ 
is the electron and n pn the phonon contribution to the 
thermal conductivity. Equation (fTJ) has a transparent 
physical meaning: The energy of the electron-hole pair 
is released as heat, see also Fig. [T] whereas the thermal 
conduction to the leads removes the heat from the grain 
until a steady state is reached. The typical energy of 
electron-hole pairs is of order of the grain temperature, 
T g , and included in the expressions for a- m and K- m . The 
leads are assumed to be a heat sink and stay at ambient 
temperature T since they are much larger than the grain. 
The first term on the r. h. s. of Eq. ([TJ) describes the 
energy released in the grain from an electron-hole pair 
times the square of the discrete gradient of voltage (V/ a) . 
The latter ensures that the result is invariant under a 
sign change of the voltage V. The second term describes 
the heat removal from the grain driven by the discrete 
laplacian of the temperature (2(T S — T)/a 2 ) with its rate 
K, cf. Eq. (23) of Rcf. lj. Our considerations are valid 
as long as the applied voltage does not break down the 
Coulomb blockade, i.e. when eV < E c . 



Before solving Eq. JT]) for temperature T g vs. time t 
exactly, we discuss its stationary solution corresponding 
to large times, ( > f , where t* is some characteristic 
time scale defined below. In other words we first esti- 
mate how pronounced the effect of inelastic cotunncling 
on grain heating is. For times t 3> t* the grain tempera- 
ture reaches some constant steady state value, T g = T*, 
i.e. dT g /dt = 0. Using Eq. {TJ one obtains 



a ia V 2 = 2 K (T* -T). 



(2) 



The inelastic cotunneling electric, a- m , and thermal, 
k, conductivities/conductance of a single metallic grain 
weakly coupled to two leads are given by the following 



expressions 
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T 2 + (eV) 2 

o 2 -1 2 g ' ^ I 
cr in = 2e a g t ^ , 

, , T 2 + (eV) 2 
Kin = l e a g t T g 
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(3a) 
(3b) 
(3c) 



In Eq. (|3b[) and (|3cp 7 e / p h arc numerical coefficients 
(7e = 32rr 3 /15 and 7ph = 8tt 2 /15, see [3), l P h the 
phonon mean free path l pn ~ a, and Qd the Debye tem- 
perature. Notice, that the temperature entering these 
expressions is the grain temperature, T g , which deter- 
mines the number of typically involved internal energy 
states of the grain 5 -. Substituting these expressions back 
into Eq. {T]) we obtain the following equation of motion 

**7^ = « [V 2 (*2 + - (3? - **) M + 7eV 2 )] , 

(4) 

where we introduced the dimcnsionlcss grain temperature 
1 g = Tg/T, time r = tS, and voltage V = eV/T, and used 
the dimcnsionlcss parameter; 
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Using these definitions, the solution of Eq. Q for di- 
mcnsionless grain temperature in the limit 1* 3> V is 
given by 



1 



(6) 



which can be simplified for small voltages V 2 <C 7 to 
1* = 1 + V 2 /j. The complete voltage dependence of 1* 
is shown in Fig. 

Time scale and estimates. — We define the typical 
timescale associated with the heating of the grain due to 
inelastic cotunncling as the time t* when the grain tem- 
perature T g reaches the intermediate value (T* + T)/2, 
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FIG. 2: (Color online) Voltage dependence of the dimension- 
less steady state grain temperature 1* g [solution of Eq. (f2J] for 
different values of the tunneling conductance g t , shown next 
to the graphs. All graphs are plotted for the typical param- 
eters given in the text and at temperature T = 10K. The 
values for V are valid for V < E c /T. 



which can be calculated in the limit 1* 3> V as 
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defined 
< 7, Eq. 



(7) 

Eq. ((6|). For small voltages, 
([7]) can be simplified to r* = 
31n2)V 2 /7) /(ja) meaning that the typical 
timcscalc t* is only weakly voltage dependent. 

Next, we will use typical experimental parameters and 
estimate the heating effect of an isolated grain weakly 
coupled to two leads with tunneling conductance g t = 

0. 1. For this we can first estimate the dimcnsionlcss pa- 
rameter 7 [Eq. (|5b[) ] for a a = lOnm metallic grain. For 
Aluminum, Nickel, or Iron the Debye temperature and 
Coulomb energy are Qd ~ 450K and E c « 1600K, re- 
spectively, which results in 7 w 7 x 10 3 . At a low tem- 
perature of T = 10K and typical voltages V R* 5mV, 

1. e. V = 5.8, we have for the dimcnsionlcss parameter 
a = 3 x 10~ 6 . Therefore the resulting increase of tem- 
perature is about IK and the resulting timescale for this 
temperature increase for a typical mean level spacing of 
IK is t* « 0.2ns (see Figs. d&E]). 

Discussions. — The typical timescale t* for the heat- 
ing process is, although rather short, still observable 
with state-of-the-art transport measurements which can 
have picosecond resolution^. However, one has to take 
into account that also sequential and elastic cotunnel- 
ing can contribute to the total electric current through 
the grain^i. Compared to the timescale t* , the typical 
life-time of an electron-hole pair (oc 1/T*) is 2-3 or- 
ders of magnitude shorter, such that our assumption of 
instant heat generation from the electron-hole pairs is 
well justified. Also the heat distribution inside the grain 
and leads due to phonons happens on a much smaller 
timcscalc than t* such that we can assume that electrons 
and phonons are in thermal equilibrium within the grain 
or lead systems. At this point we note that in general the 
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FIG. 3: (Color online) Plots of the grain temperature T g vs. 
dimensionless time t — tS. The characteristic scale Tg is 
given by Eq. ((2| and r* approximately by Eq. (0. The y-axis 
is shifted by T and scaled to 1 for better comparison. The 
graph shows plots for three different values of the tunneling 
conductance gt written next to each group of two curves which 
represent different values for ambient temperature T: T = 
10K (solid) and T = 15K (dashed). The inset shows the 
process of heating and cooling (with the same axis as the 
main plot): in the left part (r < 150) a finite voltage is applied 
[V = 5m V in Eq. Q] which is switched off at r = 150 [V = 0]. 



kinetic Eq. |T]) has an additional term proportional to the 
thermoelectric/Seebeck coefficient. However, this term is 
much smaller - at least one order of magnitude - than 
the electric and thermal conductivities, see Ref. [l8| . 

Finally, we can also study the single grain system in a 
steady state for voltages V > which is switched off to 
V = at a certain time to. In this case one should investi- 
gate Eq. (JTJ) in the absence of the first term on the r. h. s. 
with the initial condition T g (t ) = T*(V > 0) > T, where 
T is the lead's temperature. The time dependence of the 
cooling process is shown in the inset of Fig. [3l where 
the voltage V is switched off at the dimensionless time 
To = 150. The typical timescale of the cooling process is 
then given by Eq. © with V = 0, i.e. r* = In 2/(70;) 
which is longer than for the heating process. 

To conclude, we derived the kinetic equation for the 
heating of a single grain weakly coupled to two leads. 
The heating is governed by inelastic cotunneling pro- 
cesses which create electron- hole pairs in the grain. The 
heat is removed from the grain by electron and phonon 
thermal transport processes. Using this equation we ob- 
tained the steady state temperature of the grain and the 
typical timescale which is needed to reach this state. For 
a typical metallic grain of lOnm diameter at ambient tem- 
perature T = 10K the grain temperature rises about 10% 
for an applied voltage of 5mV and the related timescale 
is of order 0.2ns which is a noticeable important effect 
for low temperature thermoelectric devices like thermo- 
couples for temperature measurements. 
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